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Photon and spin dependence of the resonance lines shape in the strong coupling regime2 

Abstract. We study the quantum dynamics of a spin ensemble coupled to cavity pho- 
tons. Recently, related experimental results have been reported, showing the existence 
of the strong coupling regime in such systems. We study the eigenenergy distribution 
of the multi-spin system (following the Tavis-Cummings model) which shows a pecu- 
liar structure as a function of the number of cavity photons and of spins. We study 
how this structure causes changes in the spectrum of the admittance in the linear 
response theory, and also the frequency dependence of the excited quantities in the 
stationary state under a probing field. In particular, we investigate how the structure 
of the higher excited energy levels changes the spectrum from a double-peak structure 
(the so-called vacuum field Rabi splitting) to a single peak structure. We also point 
out that the spin dynamics in the region of the double-peak structure corresponds to 
recent experiments using cavity ringing while in region of the single peak structure, it 
corresponds to the coherent Rabi oscillation in a driving electromagnetic filed. Using 
a standard Lindblad type mechanism, we study the effect of dissipations on the line 
width and separation in the computed spectra. In particular, we study the relaxation 
of the total spin in the general case of a spin ensemble in which the total spin of the 
system is not specified. The theoretical results are correlated with experimental evi- 
dence of the strong coupling regime, achieved with a spin 1/2 ensemble. 

Keywords: Cavity QED, Strong coupling regime, Rabi oscillation, Vacuum- field Rabi 
splitting, Tavis-Cummings model, Jayenes-Cummings model 
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1. Introduction 

Quantum mechanics principles provides the fundamental basis for the development of 
future technologies which will integrate quantum processors with quantum memories. 
Electromagnetic impulses, like photons, can provide the link between the processing 
of quantum algorithms and temporary storage between computation steps. Namely, 
quantum information can be transferred by photons into a solid state memory able 
to preserve the quantum aspects of the information for a certain time. Initially, the 
coupling between atoms and resonators was first studied in atomic physics, since in 
the absence of significant atom-atom interaction, the atomic levels have large lifetimes 
compared to their solid-state counterparts. In such situation, the transfer of quanta 
between atoms and resonators can be observed. 

The coupling between atoms and a standing electromagnetic wave was introduced 
by Tavis and Cummings [TJ [2] in the 1960s, and calculated to increase with the square 
root of the number of atoms. This phenomenon has been observed for an ensemble 
of atoms [3l H] and single atoms [6] in cavity quantum electrodynamics experiments 
(cQED), in which photons and the atoms interact in a cavity. In more recent years, 
the use of a superconducting resonator with trapped molecules above it, has been 
proposed [7] to implement cQED experiments on a chip. These initial studies made 
use of the electric-field component of an electromagnetic field to coherently exchange 
photons between the field and the system under study. 

In solid state systems, magnetic two-level systems (like up and down spins) 
can have significantly longer coherence times compared to electrical analogues. It is 
therefore desirable to implement the cQED techniques developed for electrical coupling 
to magnetic coupling, although the magnetic-field component of an electromagnetic field 
offers much less coupling. This aspect is important for the implementation of quantum 
computing on a chip, using solid state materials, in which quantum information can be 
rapidly lost to the neighboring environment. Dipolar interactions between spins in such 
systems are a major cause of information loss, since it affects the spin coherence lifetime. 
But these processes can be reduced by diluting the spins in the host solid. Magnetic 
coupling between the electromagnetic field in the cavity and a spin ensemble has been 
demonstrated theoretically [8] and experimentally [HI dU QUI EG] , in the regime called 
strong coupling. This case requires a coupling strength larger than both the cavity's 
photon decay rate and the rate at which the spin losses its quantum state information. 
Aside coupling to a cavity, recent studies [131 E] show that diluted spin systems can 
be coupled with and exchange information with superconducting qubits, making them 
suitable quantum memories. 

The spin-photon coupling can be modeled by the Jaynes-Cummings (for one spin) 
or Tavis-Cummings model (for multi-spin), and it has been shown that the photon 
absorption or emission spectrum has a pair of peaks, called "vacuum-field Rabi splitting" 
|15j . We study how the energy structure changes with the number of photons in the 
cavity, and investigate how the change causes the spectrum of the admittance in the 
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linear response theory, and also the frequency dependence of the excited quantities in 
the stationary state under a probing field. 

When the number of photons is much smaller than the number of spins, the 
spectrum shows the quantum mechanical double peak structure, while it shows a single 
peak corresponding to the magnetic resonance in an driving oscillatory field when the 
number of photons is much larger than the number of spins. 

We also point out that the spin dynamics in the region of the double-peak structure 
corresponds to our experiments using cavity ringing (see also [9]) showing vacuum- 
field Rabi splitting. On the other hand, spin-photon system can show coherent Rabi 
oscillation, representing spin dynamics under a driving classical electromagnetic filed, 
in the region of the single peak structure. Between these two regions, namely when 
the number of spins and the number of photon in the cavity are similar, a chaotic spin 
dynamics appears. 

We also study line shape of the spectrum of photon absorption and emission in 
dissipative environments. We point out that the line shape is affected by several 
conditions: how many spins are involved, how strongly the system is driven, how strong 
the dissipations are, what types of dissipation exist, and whether we consider the linear 
spectrum or the spectrum for nonlinear steady state under driving force. 



2. Model 



2.1. Hamiltonian 

As system Hamiltonian we adopt the Tavis-Cummings model (a generalized Jaynes- 
Cumming model) "Hq with a driving force £ [H |2j [16] : 



u = no + Uz (i) 



where 



N / N N 



Ho = uWa + oo s St + N/2) - g I <J J] Sr + a ^ S+ 1 (2) 

i \ i i / 

Ut- = C (cte** + ae luJt ) . (3) 

Here {Si} denotes a set of two-level systems, and we call them 'spin' hereafter. Each 
is represented by the operator for a spin of S = 1/2. Here, ujq is the eigenfrequency 
of the cavity, u s is the energy gap of a two-level system, and u is the frequency of 
the driving force. Throughout the paper, we put H — 1. Here the system contains iV 
two-level systems. The dimension of the Hilbert space of the spin system is 2 N . First, 
we consider the case that only the uniform mode (i.e., the uniform sum of spins of the 
system) , 

S = f> (4) 



i=i 
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couples to the cavity, and thus the total spin of the system is conserved. (We will 
study the case in which the total spin can be changed in the section V.) Therefore, the 
Hamiltonian is block-diagonalized with the total spin. The total spin is S = N/2 when 
the number of spins is N. In the system, the sum of the photon number n and the 
magnetization, M = Yli=i Sf 

C = n + M (5) 

is conserved. Thus, the Hamiltonian is further block-diagonalized with this quantity. 
To characterize the block, we introduce a quantity n max = C + S which is the number 
of photons when all the spins are down, i.e., the state of M — —S. The size of the 
sub-block increases from 1 to N + 1 where n max is equal to 2S + 1(= N + 1). For larger 
photon numbers, the size of sub- matrix is N + 1, where extra photons exist even when 
all the spins are up (i.e., M = S). 

We adopt the basis in the form \n,M). The matrix elements are given by the 
relations 



S + \n,M) = C s , M \n,M + l), Q 



S,M 



M(M 



and 



<J\n,M) 



yn + l\n 
The matrix of T-Lq is given as 
/ 



1,M). 



(6) 
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We denote the eigenenergies in each sub-block characterized by n max as 



E(n max , k), k — 1, • • • kjj 



where 




for 
for 



< N 



^max _ 
rimax > N. 



(9) 



(10) 



We define the Rabi frequency as the energy difference (see [9] for analytical 
expression in the case n max = 1,2 and 3): 

WRabi(«max, k) = E(n m£kX , k+l)-E(n max , k), k = 1, • • • fc max (n max )-l.(ll) 

In Fig. Ufa), the eigenenergies for the case of N = 1 (i.e., the Jaynes-Cummings 
model) are plotted by • as a function of the number of sub-block characterized by n max . 
The energies are also depicted by bars on the left side. The Rabi frequencies are plotted 
by squares (multiplied by 5). There are two states, |n max , m = — |) and |n max — 1, m = |) 
in each subblock, and the difference of the two eigenstates (i.e. the Rabi frequency) is 
given by WR a bi = 2g v /n max . The bars denote the energy levels and the thin arrows 
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Figure 1. Eigen-energies (dots) and Rabi frequencies (blue squares) as a function of 
"max- (a) N = 1. The bars denote the energy levels and the thin arrows denote one 
photon absorptions (aJ). The bold arrow denotes an example of a nonlinear absorption 
((a + ) 2 ). The Rabi frequencies are multiplied by a factor of 5, for clarity, (b) N = 10: 
cigcn-encrgies and Rabi frequencies (multiplied by 50 for clarity). 

denote one photon absorption processes (a'). The bold arrow denotes an example of a 
nonlinear absorption ((a^) 2 ). 

We also give the photon number dependence of the energy and of the Rabi 
frequencies for the case N = 10 in Fig. QJb). Here, the Rabi frequencies are multiplied 
by 50 for increased clarity. Note that the Rabi frequencies have a larger distribution 
when iV ~ n. 

2.2. Photon number dependence of the Rabi oscillation 

In Fig. [2j we depict the Rabi oscillations versus time. The thin lines show the Rabi 
oscillations for the case of N = 1, i.e., the case of Jaynes-Cummings model, where 
the Rabi oscillations in each sub-block is a sinusoidal curve with a single frequency 
w — 2g v /7V ax . The bold lines show the Rabi oscillations for the case of N = 10. We 
find a simple oscillation for a small value (n max = 1), which corresponds to the vacuum- 
field Rabi splitting. We also find a simple oscillation for 77- max = 21 in which the photon 
degree of freedom behaves as a classical field, and the Rabi frequencies takes values 
close to bj = 2g\ < a > |. This oscillation corresponds to that of single spin in a driving 
field. If we see the graph more carefully, we find that the amplitude decreases slowly 
in the case of n max = 21, which is due to a distribution of the frequencies. In the case 
of riphoton = 31, the oscillation is almost perfect sinusoidal within the observed time 
because the distribution is small. For larger values of n max the oscillation looks perfect. 
Thus, we found two types of Rabi oscillations. The Rabi oscillation for small values of 
n max corresponds to the cavity ringing [3j [9], while that for large n max corresponds to the 
coherent Rabi oscillation which is found if the system is driven by an electromagnetic 

field nans]- 

Between them, at n p hoton ~ N = 10, we find chaotic behavior which is attributed 
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Figure 2. Rabi oscillations of the average number of photons, for various values of 
photon number. The thin blue curves show those for N = 1 and n max from to 21 
(starting from the bottom). The bold red curves are for N = 10 and n max = 1, 11,21 
and 31, starting from the bottom. 

to the wide distribution of Rabi frequencies. This change between the two regions can 
be used to estimate how many spins contribute to the phenomena. 

2.3. Dynamics in a dissipative environment 

To study the dynamics of the system, we adopt a simple Lindblad type master equation 
of the density matrix of the system p with damping terms [T9| 1201 [21] : 

— p = - YH, p] — k (o)ap + pa) a — 2apa^) 
dt % ' 

-l xy {S + S~p + P S + S~ - 2S-pS + ) - lz ((S z ) 2 p + p{S z f - 2S z pS z ) .(12) 

The first term describes the coherent dynamics. In the next three terms the relaxation 
effect is introduced by a conventional Lindblad form. The second term shows the cavity 
decay rate k, and the last two terms are giving the spin dephasing and relaxation, 
respectively. In general, there is an additional pumping term due to environment. 
This direct process could excite spins and thus generate emitted photons. But since 
the number of resonant bosons is generally low, we neglect these direct processes. 
Namely, we consider that the photon escapes from the system and is not created by 
the environment. We consider that the spin relaxation and dephasing is due only to 
energy relaxation processes, described here by constant values of k, 7^, and j z . 

For the Electron Spin Resonance (ESR) line shape in the case of linear response, 
we need an distribution p(0). The initial state, at the end of the pumping pulse, is 
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an excited state. In order to express the degree of excitation, we used a steady-state 
distribution function at a temperature T: 

p (0) = I e -*/*BT Z = Tle -H/k B T_ (13) 

Because the system is not in equilibrium, in this approximation the temperature T 
can be regarded as a parameter describing the degree of freedom. A more extended 
treatment of the temperature effects has been formulated in the Uchiyama's paper [22]. 
but here we take the present approximations. 

3. ESR spectrum 

The generalized susceptibility is given by the Kubo formular|19j. 

POO 

X(u) = Tr / a t [a,p](t)e-^rft = WpaM- (14) 
Jo 

Here, the last term is derived as follows [22]. Starting from the general form of the Kubo 
formula 

• poo 

X ba= lim - / dte- iuJt - et Tr[B(t),A]W, (15) 

where W is the stationary density matrix of the system, we use the following identities: 
Tr [£(*), A]W = Tr(B(t)AW - AB(t)W), and B(t) = e im B e - iHt .(16) 

Using trace properties regarding cyclic permutations, this leads to 

Tr[B(t), A]W = Tr(Be- im AWe im -Be~ iHt WAe im ) = Tr(Be~ iLt [A,W]), (17) 

where e~ lLt X is the time evolution of operator X (= e lHt Xe~ lHt ). Here we consider the 
driving force to be t;a^e~ luJt and analyze the response of the operator a. Therefore we 
have A = a, B = a' and p a [w] is the transformation of Tr e~ lLt [A, W}. With the except 
of the notation [a, p(t)](t), Eq. (j!4p is the standard Kubo formula for photon absorption. 
In the same equation, we used the following notations: 

[a,p](t) = e-* £i [a,p(0)], (18) 

where e~ lCt denotes the time evolution of the density matrix (TT2]) and 

POO 

p a [uj} = / [a,p](t)e~^dt. (19) 



o 



Here a) is a time independent operator. The notation [a, p](t) means 
e~ lHt [a, p]e~ lHt , where a and p are time independent operators. Thus can be put 
out of the integral of Eq. (TH]) . 

In the case of pure quantum dynamics, resonance frequencies of the linear response 
are given by 

w L R(nmax, k, k') = E(n max + 1, k') - E(n max , k), (20) 



Photon and spin dependence of the resonance lines shape in the strong coupling regimeQ 



where 

ft 1) " " ' j ^max(^max)] ft 1> " ' ' , ^max(^max "I - 1): 

at which the spectrum has delta-function peaks. The amplitude of each resonance is 
given by 



1[UI) — z | \'»max + J- 5 ft |« | 'imaxj ft/ | 

<H W 1 J, 



(21) 



where Z = V _ n V, e -^(«m«. 

z — '"max — *J z — 

We give the distribution of the linear response resonant frequency {wlr} for the 
case with N — 10 for each value of n max in Fig. [3J Because the interaction is described 
by a^S~ + aS + , the states with AM = 1 (and thus AE ~ uj ) can be connected directly, 
while the other resonances are realized via higher order processes. Consequently, the 
transition with wlr — has a large matrix element, while for the other resonances, the 
matrix elements are small and not visible in Fig. [3] 

In the large n max limit, we may use the semi-classical picture for the electromagnetic 
field as a — > \a\e luJt , and then the standard interaction Hamiltonian shows a rotating 
external field \a\ (S + e lult + S~e~ l0jt ). In this case, the resonance frequency is equal to 
Us, as expected for the classical ESR case. Thus Fig. [3] gives the transition from the 
quantum ringing (vacuum-field Rabi oscillation) region to the classical field region as a 
function of the photon number. 

In the dissipative environments, we adopt the quantum master equation ( ff2l) instead 
of the pure quantum dynamics. Because [a, p] (t) follows the equation of motion ( [121) 
|22j . we use its Fourier-Laplace transformation: 

iup a [u]-p a (0) = 1 [Us, PaH]-K (Z^Z Pa [u] - 2Zp a [u}Z^ + p a [u]Z*Z) , (22) 

where Z denotes a, S~ , or S z in the equation ffl~2"]) . with 

p„(0) = [a,p cq (T systcm )], (23) 

where T systcm is the temperature of the system (cavity). 

The equation (1221) is a set of linear equations of components of p a with the 
inhomogeneous term p a (0): 

p a (0) = iup a [u]+i\ [U s ,p a [u]]+K (tfZp a [u\ - 2Z Pa [u]Zi + Pa [u]tfZ) .(24) 
The explicit form of the equations for each matrix element (ij) is given by 

Pa(0)ij = iujp a [u]ij + i J i Z~2 ((^s)ifcPoMfcj - Pa[u]ik(Us)kj) 

k 

ZknPa[u] n j ^tkPa[ u ]kn^nj + P a\^\ikZ\ n Z n j 

• (25) 

\ k n J 

The ESR spectrum is obtained by the solution p a [cu] of this equation for values of 



uj as 



X(u>) = xV) - *X» = ^aVaM- (26) 
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Figure 3. (a) Resonance frequencies as a function of n max , for N = 10. The bars 
denote the strength of the matrix elements of photon absorption for the resonance 
|(i|a|j)| 2 . (b) Zoom for frequencies around luq. Only the resonances AE ~ uio have 
visible transition probability. One distinguish three regions: the quantum region 
showing Rabi splittings, the transition region and the classical region where one 
observes one single ESR peak. 



In Fig. IU we depict an example of temperature dependence of f° r N = 5, 
where the maximum value of n max is 25. The parameters g = 0.1, k = 0.01,7^ = 0.01, 
and 7 2 = 0.01 are used. There we find that has two peaks which correspond to the 
vacuum-field Rabi splitting[15] at a low temperature (T = 0.01). As the temperature 
increases, the resonances between excited states become to contribute and a single peak 
is formed. 



4. Non-equilibrium distribution of photons 

So far, we assumed that the photons in the cavity are in equilibrium at a certain 
temperature. But if we excite the cavity from outside by a finite driving force £, the 
distribution of photons can have a non-equilibrium form, as studied by Bishop et al.(16]. 
Also, nonlinear processes due to (a J( ) k ,k > 1 can contribute. The frequency of the 
nonlinear resonance is given by 

j k I;' -- ; ^max(^max) • X^J 



E(n,k)-E(n',k) 

w N L(n max , k) = 

n — n 
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Figure 4. Complex admittance for the case N = 5 and n max = 21, with g = 0.1, k 
ii.lil. 7r , y = 0.01 and 7z = 0.01. Here we plot data for T =0.1, 1, 2, 3, 4, 5, 7, 10, 20, 
50, 100, in increasing order as shown by the arrow. The line shape changes from the 
double peak to the a single peak as the temperature increases. Below the axis, we plot 
the positions of the resonant frequencies (dots); the bar below each dot denotes the 
strength of the matrix element of the resonance, as in Fig. [31 



4-1. Spectrum of a heterodyne measurement 

Under the driving force £, the system reaches a stationary state p s t which is obtained 
as a stationary solution of (1T2|) . However, the Hamiltonian fl3]) is time dependent and 
we cannot simply put dp/dt = 0. Thus, we use a rotating frame to remove the time 
dependence e ±tujt : 

U' = WHU- 1 , U = exp (iut(ata + S z )) . (28) 
This transforms quantities as 

Ua^U" 1 = ate™, US'U' 1 = S~e~ iuJt (29) 

and 

dupu- 1 _ jjdpjj-i + <EL pU -i + du P du ~ l 

dt dt dt dt 

= U^U^ 1 + iu ( (at a + S'WpU^ 1 - UpU^iata + S z )) . (30) 
dt ' 

Therefore, with this change, the equation of motion for p' = U pU~ x is given by (1T21) 

with the following form of the Hamiltonian: 

N / N N \ 

W = u a^a+u s {J2 S-+N/2)-g [atJ2 S r + a Yl S ? + ^ ( flt + a ) -^(ota+S z ). (31) 

i \ i i / 

Now, we consider the dynamics in this frame, and obtain a stationary state by 
setting dp/dt = 0: 

— [H', p st ] - n (atap st + p^ata - 2ap st at) 
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Figure 5. (a) Dependence of Q/2£ = Im(< a >)/£ on the driving frequency u, for 
different £. The blue dots show the linear response case (see Fig. 4, T=0.1) (b) Driving 
frequency ui dependence of |(a)| 2 /£ 2 of stationary states under various strengths of the 
driving force £ (c) Driving frequency ui dependence of the photon number of stationary 
states < a^a > /£ 2 under various strengths of the driving force £. The parameters are 
jV = 5 and n max = 25, with g = 0.1, k = 0.01 and j xy — j z — 0.01. The strengths of 
£ are 0.01,0.02,0.05,0.1,0.2, and 0.5, in increasing order as indicated by the arrows in 
each figure. 



-7 (S + S~ Pst + Pst S + S~ - 2S~ Pst S + )- lz ((S z ) 2 Pst + Art(£") 2 - 2S z p st S z ) = 0.(32) 

This is a homogeneous equation of p s t[IS]- The expectation values 

/ = Tr(a + a')pst, and Q = Tr(ia^ — ia)p st (33) 

correspond to the heterodyne measurement of the field quadratures. The steady-state 
transmission amplitude is given by Q. For a small strength £ = 0.001, the transmission 
shows a peak corresponding to the vacuum-field Rabi splitting. 

In Fig. EJa), we compare the line shape of the linear response x( w ) (°); an d that 
of Q/2£ for various values of £. We find that the spectrum of Q for small values of 
£ agree with x( w )- We plot the spectrum of |(a)| 2 /£ 2 in Fig. Efb). Here we find that 
x(o;) tends to move to a single peak, but |(a)| 2 keeps a double-peak structure. We also 
show the spectrum of (a)a) in Fig. [3fc). In cases with many spins and photons, there 
exist many resonant modes within the vacuum-field Rabi splitting, and the double peak 
structure is affected. 

In strong drive, the photon number shows a single peak, while the heterodyne 
signal has a double-peak shape. This can be understood from a view point of the phase 
decoherence of the photon mode. 

In the case u = ujq, 

n=0 

n=0 V "" 
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where we used the relation 

|n> = 4r( flt ) n l >- ( 35 ) 

yn! 

In the case with |'0(O)) p h = |0), the number of photons increases with time: 

a\i£t) = i£t\i£t) and (n) oc (£t) 2 . (36) 

The damping term causes a saturated value of the photon number in the stationary 
state. The dissipation term causes phase decoherence of the coherent state, due to the 
phase fluctuation 

a\i£t) = e iS (i{;t)\iZt). (37) 

By averaging over 5, the \(a)\ is strongly reduced, while on the other hand, (a^a) is 
robust against this effect. This explains the large reduction of the spectrum at its 
center u ~ Uq = Us, where the cavity photons are incoherently excited. 

5. Dissipation of the total spin 

The total spin is given by 

S(S + 1) = M 2 X + My + M 2 Z . (38) 

So far, we studied the case where the total spin is fixed to be the maximum value 
(S = N/2). This value is for the state with all the spin aligned. When the system is 
excited, the value of the magnetization M = —5+ (number of up spins) increases. In 
order to conserve the total spin, the xy components must be coherent. In the model 
Hamiltonian (j3]), the total spin is conserved and the dissipation dynamics (fT2j) does not 
affect the total spin, either. Therefore, in the present model, when the spins are excited, 
namely they are polarized longitudinally, the coherence of the transverse component 
appears. This large coherent transverse component would cause the super-radiance [23J. 

However, in practice, we may expect the decoherence of the spin state. In this case, 
the total spin is reduced. While the maximum of S 2 is given by S(S + 1) = N(N + 2)/4, 
the actual value can be reduced to the value corresponding to N phase-independent 
spins 

(<P\S 2 \<P) = (39) 

where 

l^) = 2^/2 Wu°2,---,(r N ). (40) 

{<n=±i} 

This implies that we cannot use the Tavis-Cummings model with only one value of 
S. The effective value of S may change for each frequency. Thus we must adopt a 
Hamiltonian with 2 N dimensional spin space instead of 2S + 1. 

To take into account this fact, we adopt a new type of dissipation mechanism, in 
which the total spin can change. In the process given by the master equation ffl2|) . the 
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Figure 6. Frequency dependences of (a) |(a)| 2 /£ 2 > (b) (a'a)/^ 2 , and (c)\(S 2 ) as a 
function of the frequency of the driving force with the dissipation term Eq. (|41l) for 
N = 3, ri max = 20 with g = 0.1, k = 0.01 and j xy — j z ~ 0.01 The values of £ are 0.01, 
0.03, 0.05, 0.07, 0.09 and 0.1, in increasing order as shown by an arrow in each figure. 



environment couple to the total magnetization M z = J^ i=1 Sf, and the dissipation effect 
was given by the last term in (I12p . In contrast, now we provide the thermal contact to 
each spin, and we have a new master equation: 

— p = — [H, p\ — K (a) an + pa) a — 2ana^) 
at in ' 

N N 

t=i i=i 
Using the present master equation in the extended space, we present in Fig. [6] the 
frequency-dependences of |(a)| 2 , (a^a) (normalized to ^ 2 ), and |(2M+iV)| 2 for stationary 
states. For N = 3, the maximum value of S 2 is S(S+ 1) = 3.75, and its minimum value 
is 3N/4 = 2.25 [21]. Here we find that the total spin is largely reduced, and almost 
reaches the minimum value. 

Thus, we expect that even if the number of spins is large in the cavity, practical 
values of S could be small. In the thermodynamic equilibrium state, we must estimate 
how many number of spins in the material contribute coherently (see next Section). 



6. Experimental study 

6.1. Strong coupling regime 

We have carried out a spectroscopic study of a large N = N sp i ns sample, in the strong 
coupling regime. The system of choice is the two- level spin system of the dipheriyl-picri- 
hydrazyl (DPPH), a well known Electronic Spin Resonance (ESR) standard [25]. Each 
DPPH molecule contains a single spin S — 1/2 and therefore the spins are relatively 
diluted and their environment is lacking anisotropy (g$ — 2). As a consequence, the 
line width of the DPPH resonance is very narrow, of the order of 0.1 mT. It is only 
at very high Zeeman splitting (large static fields and large frequencies, of few hundreds 
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Figure 7. Example of cavity ringing detected following a pump pulse. The readout 
is blocked during pumping, to avoid potentially large signals into the amplification 
stage. The cavity emission shows beatings which are visible as distinct peaks in the 
FFT trace (see Fig.©. 



GHz) that the anisotropy in the g-factor starts to be detectable, leading to an increased 
resonance linewidth. But at frequencies ~ 10 GHz as in our study, DPPH is among the 
best candidates to study resonance splitting consequent to cavity photons rather than 
anisotropy in spin environment. 

The sample is placed inside a cylindrical cavity operated in the TE n mode, with 
an optimized coupling to the microwave B-field. The resonance mode is located at 
Uq/2tt = 9.624 GHz and the system is operated in reflection: a microwave pulse is 
pumped into the cavity and the reflected pulse is detected by a home build heterodyne 
setup. More precisely, the reflected pulse is down-mixed by another pulse with a 
frequency close to the one used in the pump pulse (in our case, the shift is of 100 MHz 
and the pump pulse has a length of 1 (is). 

To study the eigenenergies of the coupled spin-photon system, we use the microwave 
pulse at bj to excite the system and record the photon release after the excitation. When 
the microwave is switched off, the cavity coherently emits photons corresponding to its 
own eigenmodes. The time evolution of such a signal is shown in Fig. [7] where one 
can clearly see beatings of close but distinct frequencies in the emitted signal. This 
technique is called cavity ringing and allows to study the energy spectrum by detecting 
photon frequencies in the Fast Fourier Transform (FFT) of decaying ringing oscillations. 
Obviously, most of the FFT spectrum is localized at the pump frequency. Therefore, 
we detune slightly the pump from the cavity resonance (by 50 MHz), in order to see the 
cavity-spin eigenmodes in the tail of the FFT. 

Measured Fourier traces as a function of the applied static field are presented in 
Fig. [HI The FFT intensity is gray coded, the dark areas representing the eigenmodes of 
the spin-photon coupled system. Frequencies are relative to the sample-loaded cavity 
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Figure 8. (a) Intensity plot of the Fourier transform of detected cavity ringing, 
measured at room temperature. The plot is given as a function of the magnetic field 
around the resonance condition (FFT frequencies are relative to the cavity resonance) . 
One observes a level repulsion developing in the resonance window, due to the spin- 
photon strong coupling. The effect is particularly visible at zero detuning (vertical 
marking line), (b) Actual FFT trace recorded at zero detuning showing a measured 
Rabi splitting of 10.9 MHz. The peaks are distinguishable by a splitting larger than 
peaks width. 



resonance, in absence of applied field. The static field is varied around the empty 
cavity resonance condition ij,qH z = hu>o/g s fiB, with g s = 2 by a detuning is given by 
/ioSH z = — hA/g s fi B . At zero detuning, the actual FFT plot is given in Fig. [8] and one 
clearly distinguish (that is, peaks separation is larger than peak width) two eigenenergies 
split by 10.9 MHz. The continuous lines are fits made considering a level repulsion A. 

In our experimental conditions, due to the spin-photon strong coupling, one 
observes two peaks separated by the characteristic Rabi splitting. The strength of 
the coupling can be estimated using the model of N spins coupled to a single photon, 
although the cavity contains a large number of n max photons. As we see in Fig. El as 
long as the number of photons (n max ) is much smaller than the number of spins N, 
the resonance frequencies are given by E(l,k) — E(0, 1), k — 1,2 of Eqs. <Q. Thus, 
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the difference of the resonance frequencies is given by that of the vacuum-field Rabi 
oscillation. Therefore, the cavity ringing shows coherent oscillations with the frequencies 
w e0 ,i/27r, given by: 



u e0tl -u = -^± l -^^ + n% (42) 

on which the continuous lines of Fig. Ufa) are based. The fit procedure leads to a Rabi 
splitting of Vl r /2tx = 10.9 MHz. 

6.2. Photon number and the types of Rabi oscillation 

As mentioned previously, it is important to compare the number of spins (or total 
spin size) to the number of photons, since quantum peaks require large splitting (large 
N) and a smaller number of photons. In our experiment, taking into account the 
DPPH spin density, we estimate the number of spins to be on the order of 10 20 . Since 
the temperature is T = 300K, the size of the total spin will be reduced due to the 
Boltzmann population of the two levels . For a Zeeman splitting of oj = 9.62 GHz, the 
magnetization is estimated as 

m(T = 300K ) = ^= tanh ( ~ 1(T 3 ^^, (43) 



2 \2k B T J 2 

and thus the dominant value of the total spin is of the order 

S = m(T = 300K) ~ i - 3 ^p. (44) 

The number of photons inside the cavity can be roughly estimated by using the Q factor 
of the cavity (Q ~ 1000) to get a photon life time of less than 100 ns and estimations of 
the microwave power delivered by the setup. We find a range for the number of photons 
of the order of 10 12 — 10 14 which fulfills the condition: 

"-photon < S (45) 

In this case, we expect indeed that the vacuum-field Rabi splitting of the resonance 
frequency to be clearly found. The splitting is measured to be of the order 10MHz: 

Vl k /2tx ~ gVs = 10MHz, (46) 

from which we can estimate an order of magnitude for the coupling constant g ~ 
10~ 2 - 10" 1 Hz. 

On the other hand, recent experiments demonstrate the existence of coherent Rabi 
oscillation of diluted spin systems driven by an electromagnetic field [T71 HE]- I n those 
experiments the number of photons in the cavity is much larger than N as we see in 
Fig. [21 It is therefore an interesting future problem to study the transition between 
these two regions (i.e., the region of the vacuum- field Rabi oscillation (n N) and that 
of the field driven Rabi oscillation (n ^> N). 
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7. Summary and discussion 

We studied the line shape as a function of the size of an ensemble of spins and the cavity 
mode photon number. We found the line shapes change depending on the number of 
spins, number of photon, and also types of dissipation. 

We found two regions of the photon numbers. One of them is the cavity ringing 
region where N 3> n, while the other is that of coherent Rabi oscillation under driving 
electromagnetic field. 

We presented an experimental data of the hybridization of the cavity and the spins 
which shows double peaks in the spectrum representing the vacuum-field Rabi oscillation 
due to the cavity ringing, and we estimate the number of spins and photons contributing 
to the coherent dynamics. 

With our model, the role of other factors in the spread of Rabi splitting, such as: 
dipolar or hyperfine interactions, local anisotropic crystal fields, and the size of the 
sample vs. size of a cavity (mode) can be studied as well. It is also an interesting 
problem to study the relation between this dissipation effect and the existence of the 
optical bistabilitypEJ |27l EHJ |29]. If the dissipation of the transverse coherence does 
not exist, the polarization of the system is monotonic with the driving force. But, the 
dissipation can cause a nonlinear dependence which leads to optical bistability. 



Acknowledgments 

This work was partially supported by the Mitsubishi foundation, the NSF Cooperative 
Agreement Grant No. DMR-0654118 and No. NHMFL-UCGP 5059, NSF grant 
No. DMR-0645408, National Computer Facilities (Netherlands) and also the Next 
Generation Super Computer Project, Nanoscience Program of MEXT. The computation 
in this work has been done using the facilities of the Supercomputer Center, Institute for 
Solid State Physics, University of Tokyo. We thank City of Marseille and Aix Marseille 
University for financial support. 



References 

Jaynes E. T. and Cummings F. W. 1963, Proc. IEEE 51, 89. 
Tavis M. and Cummings F. W., 1968, Phys. Rev. 170, 379. 

Kaluzny Y., Goy P., Gross M., Raimond J. M., Harochc S., 1983 Phys. Rev. Lett. 51, 1175-1178. 
Sanchez Mondragon J. J., Narozhny N. B. , and Eberly J. H. , 1983 Phys. Rev. Lett. 51, 550. 
Thompson R. J., Rempe G., Kimble H. J., 1992, Phys. Rev. Lett. 68, 1132. 

Bruno M., Schmidt-Kaler F., Maali A., Dreyer J., Hagley E., Raimond J. M., and Haroche S., 

1996 Phys. Rev. Lett. 76, 1800. 
Rabi P., et al, 2006, Phys. Rev. Lett. 97, 033003. 
Imamoglu A.,2009, Phys. Rev. Lett. 102, 083602. 

Chiorescu I., Groll N., Bertaina S., Mori T., and Miyashita S., 2010, Phys. Rev. B 82, 024413. 
Schuster D. I., Sears A. P., Ginossar E., DiCarlo L., Frunzio L., Morton J. J. L., Wu H., Briggs 
G. A. D., Buckley B. B., Awschalom D. D., and Schoelkopf R. J., 2010 Phys. Rev. Lett. 105, 
140501. 



Photon and spin dependence of the resonance lines shape in the strong coupling regimel9 

[11] Kubo Y.,. Ong F. R, Bertet P., Vion D., Jacques V., Zheng D., Dreau A., Roch J.-F., Auffeves 
A., Jelezko F., Wrachrup J., Barthe M.F., Bergonzo P., and Esteve D., 2010, Phys. Rev. Lett. 
105, 140502. 

[12] Amsiiss R. ct al., 2011, Phys. Rev. Lett. 107, 060502. 

[13] Zhu X., Saito S., Kemp A., Kakuyanagi K., Karimoto S., Nakano H., Munro W. J., Tokura Y., 
Evcritt M. S., Nemoto K., Kasu M., Mizuochi N., Semba K., 2011, Nature 478, 221. 

[14] Kubo Y., Grezes C, Dewes A., Umeda T., Isoya J., Sumiya H., Morishita N., Abe H., Onoda S., 
Ohshima T., Jacques V., Dreau A., Roch J.-F., Diniz I., Auffeves A., Vion D., Esteve D., Bertet 
P., 2011, Phys. Rev. Lett. 107, 220501. 

[15] Agarwal G. S., 1984, Phys. Rev. Lett. 53, 1732. Houdrc R., Stanley R. P., and Ilegems M., 1996, 
Phys. Rev. A53 2711. ZhuY., Gauthicr D. J., Morin S. E., Wu Q., Carmichacl H. J., and 
Mossberg T. W., 1990, Phys. Rev. Lett. 64, 2499. 

[16] Bishop L.S., Chow J. M., Koch J., Houck A. A., Devoret M. H., Thuncberg E., Girvin S. M., 
Schoelkopf R. J., 2009, Nature Physics 5, 105. 

[17] Bcrtaina S., Gambarelli S., Tkachuk A., Kurkini N., Malkin B., Stcpanov A., and Barbara B., 
2007, Nat. Nanotechnol. 2, 39. Bertaina S., Gambarelli S., Mitra T., Tsukerblat B., Muller A., 
and Barbara B., 2008, Nature (London) 453, 203. Bertaina S., Chen L., Groll N., Van Tol J., 
Dalai N. S., and Chiorescu I., 2009, Phys. Rev. Lett. 102, 050501. Bertaina S., Shim J. H., 
Gambarelli S., Malkin B. Z., and Barbara B., 2009, Phys. Rev. Lett. 103, 226402. Bertaina S., 
Groll N., Chen L., and Chiorescu I., 2011, Phys. Rev. B 84, 134433. 

[18] De Raedt H., Barbara B, Miyashita S., Michielsen K., Bertaina S, and Gambarelli S, 2012, Phys. 
Rev. B 85, 014408. 

[19] Kubo R., Toda M. and Hashitsume N., Statistical Physics II, (Springer- Verlag, New York, 1085). 

[20] Lousell W. H., Quantum Statistical Properties of Radiation, (Wiley, New York, 1973). 

[21] Walls D. F. and Milburn G. J., Quantum Optics (Springer- Verlag, 1995). 

[22] Uchiyama C, Aihara M., Saeki M. and Miyashita S., 2009, Phys. Rev. E, 80, 021128 1-16. 

[23] Dicke R. H., 1954, Phys. Rev. 93, 99. 

[24] We are preparing the massive calculation needed for a larger number of spins. 

[25] Krzystek J., Sienkiewicz A., Pardi L., Brunei L.C., 1997, J. Magn. Reson. 125(1), 207. 

[26] Lugiato L.A., in Progress in Optics Vol. XXI, ed by E. Wolf pp. 69-216. 

[27] Drummond P. D., 1981, IEEE, J. Quantum Electronics, QE-17, 301. 

[28] Gripp J., Mielke S. L., and Orozco L. A., and Carmichael H. J., 1996, Phys. Rev. A 54, R3746. 

Raizen M. G., Thompson R. J., Brecha R. J., Kimble H. J., and Carmichael H. J., 1989, Phys. 

Rev. Lett. 63 240. 
[29] Shirai T., Mori T. and Miyashita S., unpublished. 



